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We present iV-body cosmlogical simulations in the framework of the Newtonian limit of scalar- 
tensor theories of gravity. The scalar field is described by a modified Helmholtz equation with 
a source that is coupled to the standard Poisson equation of Newtonian gravity. The effective 
gravitational force is given by two contributions: the standard Newtonian potential plus a Yukawa 
potential stemming from massive scalar fields. In particular, we consider simulations of ACDM 
models and compute the density and velocity profiles of the most massive groups found at z=0. 
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I. INTRODUCTION 

X> ; 

In this work we present some preliminary results about the role scalar fields play in cosmological simulations, in 
particular on the process of large scale structure formation. Scalar fields have been around for so many years since 
the pioneering work of Jordan, Brans, and Dicke[l|, [2|. Nowadays they are considered as a mechanism for inflation^; 
the dark matter component of galaxies the quientessence field to explain dark energy in the universe 0. The 
main goal of this work is to study the large scale structure formation where the usual approach is that the evolution 
Oh 1 of the initial primordial fluctuation energy density fields evolve following Newtonian mechanics in an expanding 
background The force between particles are the standard Newtonian gravitational force. Now, we will see that 
!_i ■ we can introduce the scalar fields by adding a term in this force. This force will turn out to be of Yukawa type 
' with two parameters (a, A)Q- For so many years this kind of force, the so called fifth force, was thoroughly studied 
theoretically^ and many experiments were done to constrain the Yukawa parameters^ L We have also been studying, 
in the past years, the effects of this kind of force on some astrophysical phenomena fK ToL [Til. [l2| . The Yukawa force 
comes as a Newtonian limit of a scalar-tensor theory with the scalar field non-minimally coupled to gravitation[l3j 
although other alternatives can be found [l4j. It is our purpose to find the role these scalar fields play on the large 
scale structure formation processes. We start by discussing the standard LCDM model and the general approach in 
iV-body simulations (See Bertschinger [THj for details). Then, we present the modifications we need to do to consider 
the effects of a static scalar field and we show the results of this theory for the cosmological concordance model of a 
ACDM universe [l(| . To perform the simulations we have modified a standard serial treecode the author has developed 
[ItJ and the Gadget 1 [l8[ (see also |http : //www . astro . inin . mx/mar ) in order to take into account the contribution 
0> of the Yukawa potential. 
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II. EVOLUTION EQUATIONS FOR A ACMD UNIVERSE 



■ A. Newtonian approximation 

a 

The study of large-scale formation in the universe is greatly simplified by the fact that a limiting approximation 
of general relativity, Newtonian mechanics, applies in a region small compared to the Hubble length cH~ x (cH^ 1 s» 
3000/i -1 Mpc, where c is the speed of light, Hq = lOO/i km/s/Mpc, is Hubble's constant and h « (0.5 — 1)), and 
large compared to the Schwarzschild radii of any collapsed objects. The rest of the universe affect the region only 
through a tidal field. The length scale cHq 1 is of the order of the largest scales currently accessible in cosmological 
observations and H^ 1 w 10 10 /i _1 yr characterizes the evolutionary time scale of the universe. 

The Newtonian approximatiion can fail at much smaller R if the region includes a compact object like a neutron 
star or black hole, but one can deal with this by noting that at distances large compared to the Schwarzschild radius 
the object acts like an ordinary Newtonian point mass. It is speculated that in nuclei of galaxies there might be black 
holes as massive as 10 9 M Q , Schwarzschild radius ~ 10 14 cm. If this is an upper limit, Newtonian mechanics is a good 
approximation over a substancial range of scales, 10 14 cm <C r <C 10 28 cm. 
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B. General Scalar-tensor theory 

Let us consider a typical scalar-tensor theory given by the following Lagrangian 



L = 



16tt 



(1) 



Here t^ t „ is the metric, £m[s^v) is the matter Lagrangian and u)(<f>) and V((j>) are arbitrary functions of the scalar 
field. From Lagrangian |T]) we get the gravitational equations, 
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and the scalar field equation 
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where ()' = -S-r. The gravitational constant is now contained in V{<f>) and the scalar field get a mass tosf- 



C. Newtonian approximation of STT 



In the present study, however, we want to consider the influence of scalar fields in the limit of a static STT, and 
therefore we need to describe the theory in its Newtonian approximation, that is, where gravity and the scalar fields 
are weak (and time independent) and velocities of stars are non-relativistic. We expect to have small deviations of 
the scalar field around the background field, defined here as (<j>) and can be understood as the scalar field beyond all 
matter. If one defines the perturbations (f> = (4>) + <p and <? M „ = rj^ + h^, where rj^ is the Minkowski metric, the 
Newtonian approximation gives [l3T | 
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we have set ((f)) = (1 + o)/Gn an d a = 1/(3 + 2w). In the above expansion we have set the cosmological constant 
term equal to zero, since on galactic scales its influence should be negligible. We only consider the influence of dark 
matter due to the boson field of mass msF governed by Eq. ([5]), that is the modified Helmholtz equation. Equations 
(H]) and (O represent the Newtonian limit of STT with arbitrary potential V((f>) and function ui(<f>) that where Taylor 
expanded around (0). The resulting equations are then distinguished by the constants Gjv, a, and A = hp/msFC. 
Here hp is Planck's constant. 

Note that Eq. ^ can be cast as a Poisson equation for ip = (1/2) (/iqo + <fr/ (</>)), 



1 + a 



(6) 



The next step is to find solutions for this new Newtonian potential given a density profile, that is, to find the 
so-called potential-density pairs. General solutions to Eqs. (JSJ) and © can be found in terms of the corresponding 
Green functions, and the new Newtonian potential is 
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The first term of Eq. I|7|. given by ip, is the contribution of the usual Newtonian gravitation (without scalar fields), 
while information about the scalar field is contained in the second term, that is, arising from the influence function 
determined by the modified Helmholtz Green function, where the coupling u> (a) enters as part of a source factor. 
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D. Cosmological evolution equations using a static STT 



To simulate cosmological systems, the expansion of the universe has to be taken into account. Also, to determine 
the nature of the cosmological model we need to determine the composition of the universe, i. e., we need to give the 
values of fli for each component i, taking into in this way all forms of energy densities that exist at present. 

If a particular kind of energy density is described by an equation of state of the form p = wp, where p is the pressure 
and w is a constant, then the equation for energy conservation in an expanding background, d(pa 3 ) — — p<i(a 3 ), can 
be integrated to give p oc a~~ 3 ^ 1+w \ Then, the Friedmann equation for the expansion factor a(t) is written as 



H 2 oT^(-)- (8) 
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where W{ characterizes equation of state of specie i. The most familiar forms of energy densities are those due to 
pressureless matter with Wi = (that is, nonrelativistic matter with rest-mass-energy density pc 2 dominating over 
the kinetic-energy density pv 2 /2) and radiation with Wi — 1/3. The density parameter contributed today by visible, 
nonrelativistic, baryonic matter in the universe is $Ib ~ (0.01 — 0.2) and the density parameter that is due to radiation 
is fin ~ 2 x 10~ 5 . In this work we will consider a model with only two energy density contribution. One which is a 
pressureless and nonbaryonic dark matter with £Idm ~ 0.3 that does not couple with radiation. Other, that will be a 
cosmological constant contribution J7a ~ 0.7 with and equation of state p = —p. The above equation for a(t) becomes 
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Here, we employ a cosmological model with a static scalar field which is consistent with the Newtonian limit given 
by Eq. (J7J). Thus, the scale factor, a(t), is given by the following Friedman model, 
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where H = a/a, Q m o and £l\o are the matter and energy density evaluated at present, respectively. We notice that 
the source of the cosmic evolution is deviated by the term 1 + a when compared to the standard Friedman-Lemaitre 
model. Therefore, it is convenient to define a new density parameter by Of = + a). This new density parameter 

is such that + 57^ = 1, which implies a flat universe, and this shall be assumed in our following computations, 
where we consider (fi^, f2^) = (0.3, 0.7). For positive values of a, a flat cosmological model demands to have a factor 
(1 + a) more energetic content (f2 m and Oa) than in standard cosmology. On the other hand, for negative values of 
a one needs a factor (1 + a) less fi m and to have a flat universe. To be consistent with the CMB spectrum and 
structure formation numerical experiments, cosmological constraints must be applied on a in order for it to be within 
the range (-1,1) [3 M HE HI . 

In the Newtonian limit of STT of gravity, the Newtonian motion equation for a particle i is written as 



= -^T^ E |^ Xj .|3 - x,-|, a, A) (11) 
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where x is the comovil coordinate, and the sum includes all periodic images of particle j, and FsF(r, a, A) is 

r ^ -r/X 



f sf (r,a,A) = l + a [l + -) (12) 

which, for small distances compared to A, is Fsf(t < A, a, A) w 1 + a (l + j^j and, for long distances, is Fsf( v > 
A, a, A) ~ 1, as in Newtonian physics. 

We now analyze the general effect that the constant a has on the dynamics. The role of a in our approach is as 
follows. On one hand, to construct a flat model we have set the condition Clf n + £1^ = 1, which implies having (1 + a) 
times the energetic content of the standard ACDM model. This essentially means that we have an increment by a 
factor of (1 + a) times the amount of matter, for positive values of a, or a reduction of the same factor for negative 
values of a. Increasing or reducing this amount of matter affects the matter term on the r.h.s. of the equation of 
motion (fTTj) . but the amount affected cancels out with the term (1 + a) in the denominator of (|lip stemming from the 
new Newtonian potential. On the other hand, the factor F$f augments (diminishes) for positive (negative) values of 
a for small distances compared to A, resulting in more (less) structure formation for positive (negative) values of a 
compared to the ACDM model. For r ^> A the dynamics is essentially Newtonian. 



(c) ACDM - z=0 - SF: a=-l/2, X=5 Mpc/h ^ (d) ACDM - z=0 - SF: a=-l/4, X=5 Mpc/h 




FIG. 1: x-y snapshots at z = of a ACDM universe. See text for details. 



III. RESULTS 



In this section, we present results of the cosmological simulations of a ACDM universe with and with- 
out SF contribution. We consider the smal box initial condition in the Cosmic Data Bank web page 
(http://t8web.lanl.gov/people/heitmann/test3.html). The initial condition use a box 90 Mpc size and 256 3 
particles. It is somewhat small as a representative simulation due to lack of large-scale power but straddle a repre- 
sentative range of force and mass resolutions for state-of-the-art large scale structure simulations designed to study 
power spectra, halo mass functions, weak lensing, and so on. 

The initial linear power spectrum was generated using the fitting formula by Klypin & Holtzman[23j for the transfer 
function. This formula is a slight variation of the common BBKS fit 24j . It includes effects from baryon suppression but 
no baryonic oscillations. We use the standard Zel'dovich approximation [25j to provide the initial particle displacement 
off a uniform grid and to assign initial particle velocities. The starting redshift is z.;„ = 50 and we choose the following 
cosmology: f2 TO = 0.314 (where SIdm includes cold dark matter and baryons), f2& = 0.044, SIa = 0.686, Hq = 71 
km/s/Mpc, cts = 0.84, and n = 0.99. These values are in concordance with measurements of cosmological parameters 
by WMAP[26]. The simulations we present here use an initial condition with only 723,925 particles, obtained from 
the original initial condition by a reduction procedure based on a tree scheme. This implies that particle masses are 
in the order of 1.0 x 10 10 M Q . The individual softening length was 20 kpc//i. These choices of softening length are 
consistent with the mass resolution set by the number of particles. 

We now present results for the ACDM model previously described. Because the visible component is the smaller 
one and given our interest to test the consequences of including a SF contribution to the evolution equations, our 
model excludes gas particles, but all its mass has been added to the dark matter. We restrict the values of a to the 
interval (—1, 1) jljjl |20|. l2ll |22| and use A = 5 Mpc/h, since this scale turns out to be an intermediate scale between 
the size of the clump groups and the separation of the formed groups. 

In Fig. 1 we show x-y snapshots at redshift z = of our ACDM model. Fig. 1 (a) presents the standard case 
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FIG. 2: (a) Density profiles for one of the most massive groups at z = of a ACDM universe. The group is located approximately 
at x — 14 Mpc//i, j/ = 57.5 Mpc//i in Fig 1(a). Vertical scale is in units of p — lO lo M0/i _1 /(/i _1 kpc) 3 . (b) The corresponding 
circular velocity. 

without SF, i.e., the interaction between bodies is through the standard Newtonian potential. In (b) we show the 
case with a = 1/2, A = 5 Mpc/h. In (c) a = -1/2, A = 5 Mpc/h. In (d) a = -1/4, A = 5 Mpc/h. One notes clearly 
how the SF modifies the matter structure of the system. The most dramatic cases are (b) and (c) where we have used 
a = 1/2 and a = —1/2, respectively. Given the argument at the end of last section, in the case of (b), for r <C A, the 
effective gravitational pull has been augmented by a factor of 3/2, in contrast to case (c) where it has diminished by 
a factor of 1/2; in model (d) the pull diminishes only by a factor of 3/4. That is why one observes for r < A more 
structure formation in (b), less in (d), and lesser in model (c). The effect is then, for a growing positive a, to speed 
up the growth of perturbations, then of halos and then of clusters, whereas negative a values (a — ► —1) tend to slow 
down the growth. 

Next, we found the groups in the system using a friend-of-friend algorithm and select one of the most massive ones. 
The chosen group is located approximately at x = 14 Mpc/h, y — 57.5 Mpc/ft. The group was analyzed by obtaining 
their density profiles (Fig. 2(a)) and circular velocities (Fig. 2(b)). The more cuspy case is for a = 1/2 and the less 
cuspy is for a = —1/2. The circular velocity curves where computed using ir = G^M(r)jr. The case with a = 1/2 
corresponds to higher values of v c , since this depends on how much accumulated mass there is at a distance r and 
this is enhanced by the factor F$f for positive values of a. 

IV. CONCLUSIONS 

In general, we can say that even though we have done first numerical simulations using non-minimally coupled SF, 
the analysis we have done is insufficient to give us a clear conclusions on the role played by SF in the large-scale 
structure formation process. We will need to do a systematic study of the evolution of the two-point correlation 
function which is a mesure of galaxy clustering. We also will need to compute the mass power spectrum and velocity 
dispersions of the halos. Therefore, we will be able make sistematic comparisons with observations. However, and in 
favor of the model, the theoretical scheme we have used is compatible with local observations because we have defined 
the background field constant <<$>>= (1 + a). A direct consequence of the approach is that the amount of matter 
(energy) has to be increased for positive values of a and diminished for negative values of a with respect to the standard 
ACDM model in order to have a flat cosmological model. Quantitatively, our model demands to have 0/(1 + a) = 1 
and this changes the amount of dark matter and energy of the model for a flat cosmological model, as assumed. The 
general gravitational effect is that the interaction including the SF changes by a factor Fsf(t, a, A) ~ 1 + a (1 + J) 
for r < A in comparison with the Newtonian case. Thus, for a > the growth of structures speeds up in comparison 
with the Newtonian case. For the a < case the effect is to diminish the formation of structures. For r > A the 
dynamics is essentially Newtonian. 
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